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Abstract 


Boundary conditions involving higher order derivatives are presented for 
simulating surfaces whose reflection coefficients are known analytically, numerically or 
experimentally. Procedures for determining the coefficients of the derivatives are 
discussed, along with the effect of displacing the surface where the boundary 
conditions are applied. Provided the coefficients satisfy a duality relation, equivalent 
forms of the boundary conditions involving tangential field components are deduced, 
and these provide the natural extension to non-planar surfaces. As an illustration, the 
simulation of metal-backed uniform and three layer dielectric coatings is given. It is 
shown that fourth order conditions are capable of providing an accurate simulation for 
the uniform coating at least a quarter of a wavelength in thickness. Provided, though, 
some comprise in accuracy is acceptable, it is also shown that a third order condition 
may be sufficient for practical purposes when simulating uniform coatings. 
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Fig. 1 : Metal-backed dielectric layer. 

Fig. 2: Reflection coefficient phase vs <> for a metal-backed dielectric layer of 

thickness A/10 having e = 4 and p= 1. Comparison based on 4th, 2nd and 
1 st order boundary conditions with the "a" constants as given in (20) - (21 ): 
(a) H-polarization, (b) E-polarization. 

Fig. 3: Maximum allowed thickness vs |N| for a metal-backed dielectric layer 

modelled using the 4th order boundary conditions at y = x+, with a 2-degree 
phase (and/or 2 percent amplitude) error. Curves shown are for e = 2 and 

e = 7 with <> = 90 degrees. Results are indistinguishable for H- and 
E-polarizations. 

Fig. 4: Maximum allowed thickness vs |N| for a metal-backed dielectric layer 

modelled using the 4th order boundary conditions at y = x+, with a 2-degree 
phase (and/or 2 percent amplitude) error. Curves shown are for e = 2 and 
e = 7 with <J> .= 45 degrees: (a) H-polarization, (b) E-polarization. 

Fig. 5: Maximum allowed thickness vs |N| for a metal-backed layer modelled using 

the 1st and 2nd order GIBC with a 10-degree phase (and/or 10 percent 

amplitude) error. Curves shown are for e = 2 and e = 7 with incidence at 30 
degrees from grazing (a) H-polarization 2nd order GIBC. (b) H-polarization 
1 st order GIBC. 

Fig. 6: Maximum allowed thickness vs. |N| for a metal-backed layer modelled using 

the 1st and 2nd order GIBC with a 10-degree phase (and/or 10 percent 

amplitude) error. Curves shown are for e = 2 and e = 7 with incidence at 55 
degrees from grazing (a) H-polarization 2nd order GIBC. (b) H-polarization 
1st order GIBC. 

Fig. 7: Maximum allowed thickness vs. |N| for a metal-backed layer modelled using 

the 1st and 2nd order GIBC with a 10-degree phase (and/or 10 percent 

amplitude) error. Curves shown are for e = 2 and e = 7 with incidence at 70 
degrees from grazing (a) H-polarization 2nd order GIBC. (b) H-polarization 
1 st order GIBC. 

Fig. 8: Maximum allowed thickness vs. |N| for a metal-backed layer modelled using 

the 1st and 2nd order GIBC with a 10-degree phase (and/or 10 percent 
amplitude) error. Curves shown are for e = 2 and e = 7 with incidence at 30 
degrees from grazing (a) H-polarization 2nd order GIBC. (b) H-polarization 
1st order GIBC. 
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Fig. 9: Maximum allowed thickness vs. |N| for a metal-backed layer modelled using 

the 1st and 2nd order GIBC with a 10-degree phase (and/or 10 percent 
amplitude) error. Curves shown are for e = 2 and e = 7 with incidence at 55 
degrees from grazing (a) H-polarization 2nd order GIBC. (b) H-polarization 
1 st order GIBC. 

Fig. 1 0: Maximum allowed thickness vs. |N| for a metal-backed layer modelled using 

the 1st and 2nd order GIBC with a 10-degree phase (and/or 10 percent 

amplitude) error. Curves shown are for e = 2 and 6 = 7 with incidence at 70 
degrees from grazing (a) H-polarization 2nd order GIBC. (b) H-polarization 
1 st order GIBC. 

Fig. 1 1 : Maximum allowed thickness vs. |N| for a metal-backed layer modelled using 

3rd order GIBC with a 10-degree phase (and/or 10 percent amplitude) error. 

Curves shown are for e = 2 and e = 7 with incidence at 30 degrees from 
grazing (a) H-polarization (b) E-polarization. 

Fig. 1 2: Maximum allowed thickness vs. |N| for a metal-backed layer modelled using 

3rd order GIBC with a 10-degree phase (and/or 10 percent amplitude) error. 
Curves shown are for e = 2 and e = 7 with incidence at 55 degrees from 
grazing (a) H-polarization (b) E-polarization. 

Fig. 13: Maximum allowed thickness vs. |N| for a metal-backed layer modelled using 

3rd order GIBC with a 10-degree phase (and/or 10 percent amplitude) error. 
Curves shown are for e = 2 and e = 7 with normal incidence grazing 
(a) H-polarization (b) E-polarization. 

Fig. 1 4: Geometry of a three-layer metal-backed coating. 
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I. Introduction 


The use of non-metallic materials, possibly in the form of a non-uniform or 
multilayer coating applied to a metallic substrate, has made necessary the 
development of methods for simulating material effects in scattering. This is important 
in the analytical treatment of canonical geometries and also for the efficient generation 
of numerical solutions. 

A possible approach is to employ approximate boundary conditions [1], and 
the impedance (or Leontovich) boundary condition [2] has been widely used for this 
purpose. But this condition allows only one degree of freedom through the single 
surface impedance assumed, and there are limits to the surface properties that can be 
simulated in this manner. The inclusion of higher order derivatives of the field 
components on the surface increases the flexibility, and leads to a hierarchy of 
boundary conditions as discussed in this paper. The first order version is equivalent 
to the standard impedance condition. An example of the second order version was 
developed by Weinstein [3] and in [4,5] to simulate thin dielectric layers with and 
without a metal backing. The conditions have also been used [6,7] to simulate a 
perfectly absorbing surface in a finite element analysis of exterior scattering problems. 

In the general version considered here the boundary conditions were first 
employed by Karp and Karal [8] to study the surface waves supported by dielectric 
coatings. For a planar surface the conditions involve the normal derivatives of the 
normal field components, and it is shown here how they can be used to model the 
reflection coefficient of the surface. The required order of the condition increases with 
the complexity of the surface being modelled, and the effect of displacing the surface 
where the boundary condition is enforced is also discussed. To extend the resulting 
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conditions to a non-planar surface, they are first expressed in terms of the tangential 
components of the field, requiring the enforcement of duality. This is the same 
procedure used in the case of the standard (first order) impedance condition [1, 9], and 
it is then found that one effect of the higher order derivatives is to make the conditions 
less local in character through the inclusion of tangential field derivatives. 

To illustrate the application of these conditions, the problem of a planar 
metal-backed uniform dielectric layer is examined in detail. By expanding the known 
reflection coefficient, boundary conditions up to the fourth order are derived. It is 
shown that even for a layer thickness as large as a quarter wavelength, fourth order 
conditions applied at the upper surface of the layer provide an excellent simulation for 
all angles of incidence, and their accuracy is determined. Of course, such higher order 
conditions are not without disadvantages, and their use in analytical and numerical 
work is discussed. 

2. Boundary Conditions for a Planar Surface 

In terms of the Cartesian coordinates x, y, z, the region y < 0 is occupied by a 
laterally homogeneous material which may, however, be stratified or have properties 
varying in depth. For any field incident on this medium we seek a boundary condition 
which can be applied at the surface y = 0 and will accurately reproduce the field in 
y >0. 


I 
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The boundary conditions proposed are 



where k is the free space propagation constant and a time factor e' io)t has been 
assumed and suppressed. f m and f m are constants which are chosen to reproduce 


the desired scattering properties of the surface, and since a knowledge of E y or H y 
alone is not in general sufficient to determine an electromagnetic field, the constants 
cannot be chosen independently of one another. 

By expanding out the product factors, the boundary conditions (1 ) can be 
written in the alternative but equivalent form 



m 

m=0 (ik) 


5^0 

9y 



nM) (jk) m 


-£-h =0 

m ’ 'y 
8y 


( 2 ) 


where a m and a m ' can be expressed in terms of r m and r m ’ by equating the 
coefficients of like derivatives. As an example, for M = M’ = 3, 


a o=rir 2 r 3 , a i =r 1 r 2 +r2r3+r 3 r 1 , 
a 2 =r, + r 2 + r 3 , % = i • 
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When M = M' = 1 with H, = 1/ T 1 , the conditions are equivalent [1] to impedance 
boundary conditions for a surface with normalized impedance f”i . They can, therefore, 

be regarded as generalizations of the standard impedance boundary condition that 
allow the simulation of a greater variety of material properties through the inclusion of 
additional derivative factors. 

If the plane wave 


_i -ik (x cos<|> + y sin)>) 

E , =e 

is incident on a surface y = 0 at which the boundary conditions ( 1 ) are imposed, the 
implied reflection coefficient is 


M 

m = - n 


m =1 


r m - Sln » 

r m + sin<> 


V 1 t . x m . m 

2^H) a m sin $ 

n>0 

M 

V 1 . m 
2,a m sm * 

m «0 


(3) 


with an analogous expression for the reflection coefficient R' (<|>) associated with the 

component H y . The special case [* m = 1 , m = 1 , 2 M, constitutes a model of a 

perfectly absorbing surface [ 2 ]. We then have 


R (<))) = - tan 2M ( 4 - f] . 

and as M increases, there is an increasing range of angles about <|> = n/2 (normal 
incidence) where R is effectively zero. 
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More generally, the constants T m and T m can be determined from a knowledge 

of the actual reflection coefficient of the surface, and some possible procedures are as 
follows: 

(i) If analytical expressions for the reflection coefficients are available, 
expansion in the form (3) leads immediately to the identification of T m and T m . 

(ii) The constants can be chosen to recover the poles of the analytical 
expressions, but this may produce a less accurate simulation if the reflection 

coefficient is not a ratio of polynomials in sin<}). 

(iii) From computed or measured data for the reflection coefficients, r m and 
T m can be obtained by curve fitting. Although this could be adequate in any given 

case, it would not reveal the dependence on the material parameters of the surface. 

A few simple examples of a reflection coefficient R (<|>) are sufficient to 
illustrate these procedures. In general, a plane wave reflection coefficient can be 
written as 


m 

T|(<j)) + sin<)> 


( 4 ) 


where q (<()) is an angle-dependent surface impedance. If r| = rj 0 , independent of 

angle, the corresponding boundary condition has T 1 = q 0 with T m = 0 for m > 1 , and this 
is the standard impedance boundary condition [2]. Alternatively, if 
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then 


Ti(<t>)=n n + n 


• 2 X 

sin <|> 


with T m = 0 for m > 2. As another example of procedure (i), if 

— cos<|) , 

then for <j> < tc/4 we can write ti(<|>) = 1-1/2 sin 2 <> without significant loss of accuracy, in 

which case T 1 = 0.732 and \ 2 = -2.732 with f m = 0 for m > 2. For <\> > rc/4 we can set 

a = it/2 - 0 and use the same approximation for cosa. Alternatively, using procedure 

(ii), it is found that if 0 < Re. <|> < 7t/2, R(<j>) has a single pole at sin<> = 0.707 = - l"i . The 
resulting first order boundary condition is obviously not as accurate as the second 
order one given by procedure (i). 

In addition to modelling a reflection coefficient, the generalized boundary 
conditions (1 ) allow some flexibility in the location of the surface where they are 

enforced. To show this, consider a surface y = 0 at which the reflection coefficient is 

R(<j>). The corresponding reflection coefficient at y = x is then R(<|>) exp (2ikx sin<», and 
since 
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1 2 2 

1 + ikx sin<(> + 7 T (kx) sin <)>+.... 

2iktsin^ 2 

e = - - - , (5) 

1 - ikxsin<|> + — (kx) sin <}>-.... 

the phase factor can be simulated using additional derivative factors in the boundary 
condition. In particular, if kx is so small that terms O {(kx) 2 } can be neglected, one more 

derivative factor suffices with [ 1 = i/kx, whereas to the next order in kx two derivatives 

are needed with T 1 2 = (1± i)/(ikx). This allows us to consider separately the modelling 

of the reflection coefficient and the phase factor, but if the location of the simulating 
surface is not specified a priori, it can be chosen to minimize the order of the boundary 
condition for a given accuracy. 

3. Equivalent Forms of the Conditions 

The boundary conditions (1) are scalar conditions in as much as each 

involves only a single field component, but in the first order case when M = M' = 1 with 
17 = 1/T 1 they can be expressed [9] in terms of the tangential field components as 

E„.r,ZH 2 . Ez--r,ZH„ 

implying (6) 

A A A 

n x (n xE ) = - 1, Z n xH 

A 

where n is the unit vector outward normal to the surface. This provides the needed 

extension to a non-planar surface and is the form in which an impedance boundary 
condition is usually stated. In a similar manner, any pair of generalized boundary 

conditions can be expressed in terms of the tangential field components provided T m 
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and f m ' are appropriately related, and to illustrate the procedure, we consider the 


second order case. 

The boundary condition for E y is 

. -^•+ik(r 1 +r 2 )^--k 2 r 1 r 2 \ E = 0 , 

lay ay j y 

and using Maxwell's equations and the fact that V • E = 0, this can be written as 


JL 

dx 



I~1 l~2 + 1 

fi + r 2 


ZH 


1 

+ iMf, +r 2 ) 



= - -2- 


9z 


^2 + 1 i a^, 

z r 1 + r 2 ik^+fj) a z 



for any function f = f (x,z). Similarly, from the boundary condition for H y , 


(7) 


_9_ 

9x 




1 

ik ( r; + r 2 ) 



, . r; r 2+ 1 

= . -d-lH - . -- Y E.. + 


9H. 


9z 


fi + r 2 


ik(r; + r 2 ) az 


ag 1 

dx 


and therefore 


_9_ 

9z 



r;-r 2 

fi f 2 + 1 


ZH z + 


+ [2 9g 

r;r 2+ i ax 
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2 2H„ + -!^i-Z^ + 


r, r 2 + 1 & sk ( r 1 r 2 + 1 ) 



(8) 


for any function g = g (x,z). Choose now 


f=- 



and 


l"i I" 2 + 1 1 

g = — YE . 

r 1+ r 2 iMr 1+ r 2 ) y 

Then, if 


^1 ^ [2 + ] 

l"l l" 2 + 1 ^1 + ^2 

(7) and (8) imply 


0) 



jL 

dz 



r /2 + 1 1 

— — ZH z + , 

r 1+ r 2 *k (f 1 + f 2 ) 




-9_ + JL_ 
tax 2 az 2 j 



Ti f 2 + 1 ■] 

■ 7 r -ZH x + r • 

r 1+ r 2 ik(T 1+ r 2 ) 


aE y ! aH v 

— + rA z- 1 

ik^Tg + l) dx 


= 0 , 


and as shown in [2], the only allowed solutions of these equations are zero. 
Hence 
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( 10 ) 


r, r 2+ i 

E = - .. ZH.- 


8E., 


3K 


r 1 +r 2 2 iMr 1 + r 2 ) ax ik(fir 2 +i) & 


Ti r 2 + i 

E =- ; 2 . ZH • 


1 *Ey 1 ^H y 


r 1 + r 2 x iMi^+fg) dz ik (f., r 2 + 1 ) & 


(id 


on y = 0. We note that the relation (9) connecting the f m and f m is equivalent to 


+ a Q 


a 2 + a o 


and is a consequence of duality. 

The vector form of (10) and (11) is 


( 12 ) 


n x 



E] 


a 2 + a 0 


A 

Zn x 


H + ij ^V(n.Hl 


( 13 ) 


and this provides the extension of the second order conditions to a non-planar surface 
since (1 3) is now geometry independent. Similarly, in the third order case we obtain 


A 

n x 


nxU+ „ ■ 1 ■ V 


ik(a 3 +a 1 ) 


ajn-E-^VE 


a 2+ a 0 

a 3 + a l 


Znx\H + 


ik(a 3 ' + a^) 


a ? ' n • H - tJ- V • H 
1 ik s 


(14) 


provided 
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a 2 + a 0 


( 15 ) 


83' + a/ 

+ a 0 a 3 + a i 

or, equivalently, 

I I « III 

i"i +r 2 +r 3 +r 1 r 2 r 3 
Ti r 2 + r 2 r 3 + r 3 ri + 1 

where V s is the surface divergence. The extension of (13) and (14) to still higher order 

boundary conditions is obvious. It is evident from the above that one effect of going to 
higher order boundary conditions is to make them less local through the inclusion of 
tangential derivatives of the fields. 

4. Metal-Backed Layer 

A geometry of practical interest is a uniform dielectric layer of thickness x 
backed by a metal (see Figure 1 ), and we now seek a simulation using generalized 
boundary conditions of the form (1 ) or (2). In doing so, particular attention is given to 
layers of reasonable thickness since, as discussed below, approximate boundary 
conditions are already available for very thin coatings. For this geometry the exact 
reflection coefficients are known, and we can therefore derive the conditions using 
procedure (i). These can be referred either to the surface y = 0+ of the dielectric layer 

or, by using the expansion (5), to the surface y = -x+ of the metal backing (note that (5) 

is applicable here only after letting x -x). 

1 1 


i~] r 2 + r 2 r 3 + r 3 i~i + 1 
(*i +r 2 + r 3 +r 1 r 2 r 3 


For an H-polarized plane wave having E y * 0 the reflection coefficient is 


m = - 


7 N 2 - cos 2 (|) tan - cos 2 (|) ) - ie sin<p 
7 N 2 - cos 2 <t> tan (kt J N 2 - cos 2 <|) ) + ie sin<t> 


(17) 


where N = V ep is the complex refractive index, e is the relative permittivity and |i 

is the relative permeability of the layer. The corresponding reflection coefficient for 
E-polarization is 


R’ ((D) = 


7 N 2 - cos 2 <p cot (ktyV - cos 2 <p ) + ip. sin<p 
7 N 2 - cos 2 <D cot (kt 7 N 2 - cos 2 <|) ) - ip sin<|> 


(18) 


where R(<D) and R' (<D) are both referred to y = 0+. Expansion of the tangent and 

cotangent in powers of sin<D leads immediately to the identification of the constants 
a m and a m , but before doing this, it is of interest to examine two boundary conditions 

already available in the literature. 

The most commonly-used boundary condition for simulating a metal-backed 
layer is the standard impedance condition [1]. This is a first order one which can be 

I 2 2 ~ 

derived from (15) and (16) by writing yN - cos <D = N, giving 


T 1 = 7 = - i —tan (Nkx) 

fi e 


(19) 
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with all other T m = 0. For kx < 0.52 the conditions can be transferred to the surface 

y = -x+ with a maximum phase error of 10 degrees by employing the first two terms in 

(5). The resulting boundary conditions are then second order ones in which T 1 and 
T 1 remain the same and \ 2 and \ 2 = (ikx)’ 1 . These have been used [4] to simulate a 

metal-backed dielectric half-plane and since they are valid only if |N| » 1 , they will be 
referred to as high contrast conditions. Low contrast conditions can be derived by 

introducing the approximations tan x » x and cot x * 1/x in (17) and (18), giving 


RM>)- 


2 2 
kx sin <J> - ie sintf) + kx (N - 1 ) 

2 2 
kx sin <|) + ie sin<t> + kx (N - 1 ) 


and 

* w . "" I *** 1 , 

ikx p sin<j> - 1 

from which the constants a m and a m are easily found. The corresponding boundary 

conditions can be transferred to the surface y = -x+ using the first two terms in the 

expansion (5), and when only the terms of leading order in kx are retained, the 
constants are 


a 0 = kt(N 2 -1) 8,-1 

I 

a 1 = ie a 1 = -ikx (p- 1) 


a 2 = -kx (e - 1 ) 


a 2 = 


0 . 


( 20 ) 
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These satisfy (11) and the resulting boundary conditions are identical to those derived 
by Weinstein [3] by expanding the fields in the dielectric as Taylor series in y. 

The high and low contrast boundary conditions provide an accurate 
simulation of the coating over the entire angular range only for limited ranges of |N| 

and kx. However, by going to a higher order condition, it is possible to produce a 

simulation that is valid for all |N| and a wider range of kx. To this end we write 

F7 , 2 2~ kl 1 sin 2 <{> 

VN * cos <]> =N- • 

By also employing the approximation 
I kx . 21 kx . 2 X 

,an l2N Sln, fj*2N Sln * ’ 

R(<|>) can be expressed as 

4 

2,H) a m sin $ 

RW-- 2 ^ (21) 

X a m si M’ 

m-0 


where 



a 2 = 2 ^{ tan ( ktN ) ' tan (^r) + kx (N - ^j-) [l + tan (kxN) tan (■—)] 


2N 


2N 


2N' 


( 22 ) 


ikxe 


a 3~' 2N 


tan (kxN)-tan (|^) 


a 4 = 


kx_ 

4N 2 


1 + tan (kxN) tan (~) 


Similarly, for E-polarization, 


a o = (2N ' 1 ) 


a 1 = i 2N|x 


1 + cot (kxN) cot (^j-) 


kx 


cot (kxN) - cot (— ) 


kx 


a? = 1 + cot (kxN) cot (—■) + kx (N - 


a 3 = -i kxp 


• __kx_ 
a 4" 2N 


1 + cot (kxN) cot (^) 


kx 


2N 


cot (kxN) - cot (^) 


cot (kxN) - cot (— ) 


(23) 


The boundary conditions implied by (22) and (23) are referred to the surface 
y = 0+ and are fourth order ones which satisfy duality. As expected, their accuracy 


improves with increasing |N|/(kx) and in Figure 2 it is shown that they predict the 


correct reflection coefficient for e = 4, ji = 1 and x = 0.1X, i.e., Nkx = 1 .26. When we set 
a 3 = a 4 = 0 and a 3 = a 4 = 0 the resulting boundary conditions are second order ones. 

As evident from Figure 2, their accuracy is substantially less, and the standard 
impedance (first order) conditions are valid only for incidence close to grazing or 
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normal. 



The accuracy of the fourth order boundary conditions is quite remarkable 
even for fairly thick coatings. Since the conditions correctly predict the dominant 
surface wave modes (see comment below), their accuracy is greatest near grazing 
incidence, but even at other angles the phase error is less than 2 degrees with 

coatings up to X74 thick for either polarization and material composition provided 

|N| > 2. This is illustrated in Figures 3 and 4 where, for e = 2 and 7 with p real, the 
maximum layer thickness for a 2-degree phase error is shown as a function of |N| for 

<f> = 90 and 45 degrees, respectively. It should be noted, in contrast to the standard 
impedance boundary condition, the accuracy of the conditions specified by the 
constants (22) - (23) is least at normal incidence and increases with the refractive 
index and loss in the material coating. 

If a 10 degree error in phase and/or 10% in magnitude is acceptable, a truncated 
form of the fourth order conditions implied by (22) - (23) become useful. Figures 5 to 
1 0 show the maximum thickness for which the second order (truncated from the fourth) 
condition is capable of predicting the coating's plane wave reflection coefficient within 
10 degrees of its actual phase and/or 10% of its actual magnitude. As seen, in 
comparison with the standard impedance boundary condition, the second order 
condition provides substantially better accuracy for incidence angles away from 
normal. Notably, the simulation improves monotonically as one approaches grazing. 
We may conclude from Figures 5 to 7 that for H-polarization, the second order 
condition is capable of simulating coatings having thickness up to 1/5 of a wavelength 
for incidence angles greater than 35° from normal (55° from grazing). This is 
regardless of the dielectric's properties since the simulation improves substantially as 
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N and/or the loss in the coating increases. In contrast, the 1 st order (standard 
impedance) condition provides a superior simulation with respect to the second order 
condition for the rest of the angular region (i.e. within 35° from normal). Turning to 
Figures 8 to 10 one again arrives at similar conclusions for E-polarization. However, it 
should be noted that for small N the deterioration of the simulation provided by the 
second order condition as normal incidence is approached, is now more rapid. 

Figures 11 to 13 show the maximum thickness for which the third order condition 
is capable of predicting the coating's plane wave reflection coefficient within 10 
degrees and/or 10% of its magnitude. It is observed that this condition provides an 

acceptable simulation for coating thickness of at least 0.4X regardless of material 
properties, angle of incidence and polarization. 

For the standard impedance boundary condition, results analogous to those in 
Figure 4 are shown in Figure 14. Even with the allowed phase error increased to 10 
degrees, the maximum layer thickness is substantially less except at normal incidence, 
particularly if |N| is small. In addition, the accuracy of this boundary condition 
decreases as grazing incidence is approached. 

Higher order boundary conditions for a metal-backed layer can also be 
derived using procedure (ii) of Section 2, which requires a knowledge of the complex 
poles of the reflection coefficients (17) and (18). The subset lying in the proper half of 
the complex plane are the usual surface wave poles, and the implied expansions of 


the reflection coefficients are 



si 

RM>) = -2,- 


sin 0 m - siruj) 


m-o sin 6 m + sin<j> 


M* 

R (<W - • S 


sin 0 m - sin<J> 


m-o sin 0 m + sin<)) 


(24) 


where f = -ik sin 0 m and f = -ik sin 0 m are the propagation constants for the 

m mm m r r 9 

H-polarized (TM) and E-polarized (TE) surface waves, respectively, provided 

Re. ( f m- f m) > 0 • 

It has been observed that the accuracy of the approximate boundary 

conditions deteriorates significantly as <}> decreases unless the correct dominant 
surface wave fields are predicted, and the boundary conditions implied by (22) and 

(23) must therefore correspond to the correct surface wave poles. The number of 
poles (or zeros) of the reflection coefficients (1 7) and (1 8) depends on the value of 

Re. (kxN). In the TM case with kxN small, only one pole exists corresponding to the 

lowest order surface wave mode. As kxN increases, two additional poles appear, one 

of which is associated with a surface wave mode when u > tc where 

u 2 = (kx) 2 (N 2 - 1 ) - sin 2 <j>. In the TE case, no poles exist for small kxN, but as kxN 
increases, two poles appear simultaneously, and one of these is associated with the 

lowest order TE surface wave mode when u > tu/ 2. An examination of the accuracy of 

(24) showed that they provide a good simulation of the reflection coefficients for those 
values of kxN such that the corresponding poles (with sin 0 m , sin 0 m not much greater 


I 
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than unity) have just appeared; however, the accuracy decreases substantially as kxN 
increases beyond the mid point between the last included pole and the next one to 
appear. 


5. Multilayer Metal-Backed Coating 

Of more practical interest is the simulation of multilayer coatings and in this 
section we consider the derivation of higher order boundary conditions applicable to a 
three layer coating with arbitrary constitutive parameters. The geometry of the coating 
is illustrated in Figure 15. For E-polarization (H y * 0) the plane wave reflection 
coefficient referenced to y = 0 + can be written as 


R' (♦) = 


f n W 

f d (<d) 


(25) 


where 


F n (d>) = -ip^Hs k y0 kyiky 2 sin (k^) cos (k^x.,) cos (k y1 x t ) 

D 

+ kyo4 ^ (k y3 h) sin (ky2 X 2> sin ( k y1 X l) 

-'>#2 C0S (ky3 X 3> S ' n (ky2 X 2) C0S (kyl X l) 

-'\h\h kyO ky 2 ky3 C0S (ky3 ^ 3 ) C0S (ky 2 X 2> Sin (kyl * 1 ) 
2 

±M 2 M 3 kyl ky 2 S ' n (ky3 * 3 ) C0S (ky 2 X 2 > S ' n (^1 X l) 

±Hl M 3 Kyl 4 S ' n (ky3 X 3) S ' n (ky 2 X 2) C0S (kyl 
2 2 

±M 2 kyj ky 3 cos (ky 3 x 3 ) sin {k# x 2 ) sin (k^ x,) 
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'W 2 kyl ky2 ky 3 C0S 0^3 x 3> C0S (ky2 X 2> C0S 0^1 X 1> 


(26) 


In the above (p. n , e n ) are the relative constitutive parameters of the nth layer with n = 1 
corresponding to the top layer, t n denotes the thickness of the nth layer, and 

kyn = k V N n ~ COs2 ^ (27) 

with N n = E n and k^ = k sin<|>. As before, <|> is measured from the surface of the 

coating. 

Similarly, for H-polarization (E y * 0) the corresponding reflection coefficient 
referenced to y = 0+ is given as 

F N (<t>) 

R (<(,) = — — (28) 

f d (« 

where now 

F n (cJ)) = k^ k y1 ky 2 cos (k^ x 3 ) cos (ky 2 x 2 ) cos (k y1 x^ 

D 

2 2 

-£^3 kyo ky2 cos (ky3 x 3 ) sin (ky 2 x 2 ) sin (ky, xj 

* e 1 e 2 kyo kyl ^3 sin (ky 3 x 3> sin (ky2 x 2 ) cos (kyi x l) 

* e l e 2 kyO ky2 ky3 Sin (><y3 X 3> C0S ^2 X 2> S ' n (kyl x l) 

±ie 2 e 3 kyi ky 2 cos (ky 3 x 3 ) cos (k^ x 2 ) sin (k y1 x,) 

2 

±,e i e 3 kyi ky 2 cos (ky 3 x 3 ) sin (ky 2 x 2 ) cos (ky, x,) 

+' 4 *Vi ky3 sin (ky 3 x 3 ) sin (k^ x 2 ) sin (k y1 x^ 

±ie i e 2 kyi ky 2 ky3 Sin (ky3 * 3 ) C0S (ky2 X 2 ) C0S (kyi (29) 
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To find the constants a m and a m appearing in (2) and (3) we must expand the 
terms sin (k yn x n ) and cos (k yn x n ) in powers of sin<j>. The simplest case is to assume 

that k yn t n is sufficiently small so that we may set 


sin 

(30) 

and 


cos (ky n T n )*1 

(31) 

By retaining only terms of 0 (x n ) we obtain 


(N^-l)kx, (N|- 1 )kx 3 

Sq - + + 

01 £2 £3 


a 1 = +i 

(32) 

. (h x 2 h \ 
a 3 = k V + - +-/ 
e i £2 £3 


and 


ao = 1 


a, = -ik (p 1 t 1 + \s^x 2 + 113 X 3 ) 

(33) 

83 = 0 . 


The last set of constants imply a first order condition, but by retaining all terms of the 


expansion we find that 

• 2 2 
ao = P1P2 ' a i (^1 " 1 ) * °2 (^2 ""0 

a, = -ik (h x i + W2 + ^3) ( 34 ) 

a2 = - {a, + a,) 

in which 
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and 


(35) 


a i = k 2 M-1M3 Ph + 

2 

a 2 = k M2 x, [ x 2 + 103 X3] . 

It is obviously expected that (34) will allow the simulation of thicker layers than (33). 

By their derivation, the boundary conditions implied by the constants in (32), 
(33) or (34) are applicable to very thin coatings ( k ^ t n < 0.6) with moderate values 

of N n . However, as noted earlier for the single layer coating, higher order conditions 
are required to allow simulations which remain valid for a wider range of N n and kt n . 

To derive a higher order condition for the three layer coating, we may parallel the 
approach employed for the single layer coating and set 

k^k^-cos 2 * -k(N n - 2^ + fj^-) (36) 

This, however, will lead to 12th order condition that is obviously impractical to employ 
analytically or numerically. Instead, a more reasonable approach is to assume that the 
three layers comprising the coating have varying refractive indices. In practice the top 
layer is low contrast dielectric, whereas the bottom layer may be a high contrast 
dielectric. With this assumption we may then set 

sin (ky, 

cos (k y1 Ti) * 1 

K y2 = k7N 2 -cos 2 *-k(N 2 -^- + |l j ^) 
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-i i 2 (IV 1) N 3 x 1 tan [(N 2 - 2 ^-) k x 2 ] tan (N 3 k Xg) 
■ e 1 h. ( N 2 ‘ N 3 tan ( N 3 k 

= - e 1 eg e 3 (n 2 - ) 

+ e? e 3 (Ng - 1) tan [(n 2 - ) k xj>] 

+ 4 N 3 tan [( N 2 ' 2^) k ^2] tan ( N 3 k x 3) 

+ £2 (nVK x 1 tan (N 3 k x 3 ) 


- *i £2 e 3 (^1 " "*) (^2 ' 2N 2 ^ Ti 2N 2 ^ " 2N 2 ^ ^ Tg 


u 


+i £, 


1 2 1 
(^2 * 2N 2 ^ X 1 + ’ e 2 e 3 (^1 2N 2 ^ 


2N 2 
k x. 


+ie n £3 (n|- 1) 2N^i£ 1 £ 3 tan[(N 2 -— ^)k xj 

K T 

-i£ 2 (N* -1)N 3 Ti 2 j^an (N 3 k x 3 ) 

-i £ 2 N 3 x, tan [(n 2 - k x 2 ] tan (N 3 k x 3 ) 

e, £2 ( N 2 ' 2^) N 3 ^ tan [( N 2 - 2^) k J tan ( N 3 k x 3 ) 

+i e 1 £2 • N 3 tan (N 3 k Xj) 


a 3 -E , c 2 e 3 2 ^ + e 1 e 2 E 3 (N 2 - I ^)^-ta n [( N 2 - i ^)kT 2 ] 


+ e i H h »" [(n 2 - ^p) k tj + 4% (Nj - 1) 


k t 2 
2 Nl 


2 . . k *2 


+ £ 1 ®2 N 3 2FT tan < N 3 k X 3> 


* e i h (^2 ‘ 2N 2 ^ Tl 2N 2 * an t(^2 ' 2N 0 ^ ^ 2 ] * an (^3 k * 3 ) 


+ e i e 2 2 KT N 3 x i tan(N 3 k h ) 


(43) 


a 4 = 4 e 2 e 3 (N^ 1 )^x 1 ^tan[(N 2 — )kt 2 ] 


*' e 2 h ( N 2 - 2 n ) Tl 2N 0 ten ' 2N, ^ k 


1 k X 2 
■h e 9 £o t< + 1 e, £o -T77- 

2 3 2 n 1 1 3 2N P 


k x 

•' e 2 N 3 X 1 2N~^ an ( N 3 k x 3> 


' £l ^ 2N ^ 3 2N ^ t ' (^2 ‘ 2N ^ ^ ^ 2 ] * an (^3 k ^ 3 ) 


(44) 




2N 2 2N 2 


k T 0 


+ e i^> t i IFE 
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(45) 


' e 1 e 2?RT N 3 t 1 ?f^ tan [( N 2 ' oRT") k tan ( N 3 k X 3> 


2No ” 3 ’ 1 2N,““"i v ’ z 2N, 




2N 2 1 2N 2 


Similarly, for E-polarization, we find 


(46) 


a’o = ^2 ^3 ( N ? ■ !) ( N 2 • 2 ^) x i tan ( N 3 k x 3) 

+ m ^3 (N 2 - 1 ) tan [(n 2 - k x 2 ] tan (N 3 k x 3 ) 
+ m|(N?-1)N 3 x 1 tan [(n 2 .~) k tj 

•h ^2 ( N 2’2^) N 3 


(47) 


a; = + i m n 2 Ji 3 (n 2 - 2 ^-) tan (N 3 k x 3 ) 


- i h M 3 ( N 2 ’ 1) x i tan [(N 2 - k 0 x 2 ] tan (N 3 ^ x 3 ) 
+ i M-i M -2 N 3 tan [(N 2 ~)koX 2 ] 


+i h ^ ( n 2 - 2^r) n 3 x i 


(48) 


' 2 1 
32*M2M3( N 1 - 1 ) 2N~ X i tan ( N 3 kx 3) 


+ ^ ^ ( N 2 ’ 2fr) T i tan (NsK,^ 
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-H 2 \h (Ni 2 - 1) (n 2 - x, 2 j^ tan [(n 2 - --) kx 2 ] tan (N 3 k x 3 ) 

• + m tig tan [(N 2 - ^-) k x 2 ] tan (N 3 k x 3 ) 
k t 

+ Mi 43(N 2 -1) 2j^-tan (N 3 kx3) 

+ A N 3 x, tan [(N 2 - gjq-) k x 2 ] + m| (n| - 1) N 3 x, 

' h M 2 ^ N 3 + h M2 (N 2 - 2^) N 3 tan [( N 2 - 2 ^) k h] (49) 

a 3 = - i m ^ 2 ng (N 2 - ^ tan [(n 2 - ^-) K x 2 ] tan (N 3 k x 3 ) 

2^7 ten ( N 3 kT 3) 

2 k T 

-i \i] Mg (N 2 - 1) tan (N 3 k x 3 ) 

-i |ii ^3 x, tan Ite-s; ") k x 2 ] tan (N 3 k x 3 ) 

2 kx 2 
^hH2N 3 — 

- j A M2 ( N 2 - 2 ^) N 3 x i ^-tan [(n 2 - ^-) k x 2 ] 

+i Mi 2 M 2 2^ N 3 Xi (50) 
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a 4 = - Ma (N? - 1) 2 ^-x, 2 ^- tan [(n 2 - ^p) k x 2 ] tan (N 3 k 


h) 


• M 2 \h ( N 2 * 2 Fr) x i tan [(N 2 - ^Jpk x 2 ] tan (N 3 k 


h) 


1 kx 2 
+ ^ ^ 2N~ Xl tan (N3 k Ts) + ^ ^ 2N~ tan (N3 k x 3> 


2 

K, x : 


2 * S )^2 1 k x 2 T / 1\1 

^ N3 x, — + M, — n 3 — 'a" l(N 2 - gjj-) k t 2 J 


(51) 


-i h 1*2 Megjp tan ~~ tan [(N 2 -^p) kx 2 ] tan (Ngk^) 


2 k T 

'' f'l ^3 X 1 2^“ tan ( N 3 kx 3) 


-i ^ ^ is: N 3 1 < ,an [ < n 2 • is:) k J 


(52) 


% = * ^2 Mg 2 ^ x i tan [ (N 2 - k t 2 ] tan (N 3 k tg) (53) 

A detailed examination of the accuracy of the generalized boundary conditions implied 
by the constants (40)*(53) as a function of thickness and constitutive parameters is, of 
course, impractical because of their dependence on numerous geometrical and 
material parameters. However, guidelines on their accuracy can be drawn from the 
single layer simulation presented in the previous section. 
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6. Concluding Remarks 

As illustrated in the case of a metal-backed uniform dielectric layer, the 
generalized boundary conditions (1) provide an excellent simulation of the scattering 

properties of the structure for thicknesses up to X/4 or more for all material properties, 
directions of incidence and polarizations. For any other structure whose reflection 
coefficient is known analytically, the appropriate boundary conditions can be derived 
in a similar manner and their accuracy quantified. Alternatively, if only computed or 
measured data for the reflection coefficient are available, the required boundary 
conditions can be found by curve fitting. 

For a boundary condition of any given order, the accuracy achieved 
depends on the location of the surface where the boundary condition is applied, and 
one advantage of the method we have followed is that the location is treated 
separately. In the example discussed, the simulating sheet was placed at the upper 
surface of the dielectric, and a fourth order boundary condition was found to produce 
excellent results. For other locations, it is necessary to expand the phase factor as 
indicated in (5), leading to additional derivative factors in the boundary conditions. 

The optimum location minimizes the complexity of the boundary conditions and/or the 
error in simulation, and the choice of another location could limit the accuracy 
achievable. If, in the above application, the simulating sheet was placed at the metal 

backing, the maximum layer thickness that can be reasonably handled is of order 0.1X. 
With the sheet located in this manner, the result of coating the metal is simply to 
replace the perfectly conducting boundary condition with the appropriate generalized 
one. It is then only a trivial extension to simulate a metal coated on both sides. 
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Provided the constants l* m and f m (or a m and a m ) satisfy a duality relation 

(see, for example, (15) or (16)) the boundary condition can be expressed in terms of 
the tangential field components. This shows that one effect of the higher order 
derivatives in (1 ) is to make the boundary conditions less local in character, and the 
resulting "vector" conditions provide the natural extension of (1 ) to non-planar 
surfaces. 

In spite of the apparent complication of the generalized boundary conditions 
and their vector equivalents, it has been found possible to work with them numerically 
and analytically. Regardless of the order of the conditions, a sheet subject to them 
supports only tangential electric and magnetic currents. Thus, in a numerical solution 
of a scattering problem, the number of unknowns is the same as for the standard 
impedance (first order) boundary condition, and to better simulate coated surfaces, we 
have already begun to incorporate the higher order conditions into our existing sheet 
scattering codes employing either the moment or conjugate gradient FFT methods 
[5,10]. 

Analytically, it is important to be able to determine the diffraction coefficient 
for the edge of a half plane or wedge subject to these boundary conditions, and 
thereby extend the capability of GTD scattering codes. In the case of a half plane we 
can use either the Maliuzhinets or Wiener-Hopf techniques and [4] is an example of 
the application of the former, whereas in [11,12] the Wiener-Hopf method was 
employed. Provided care is taken to ensure that all Fourier transforms exist in the 
classical sense, the Wiener-Hopf method can handle generalized boundary conditions 
of any order [1 3]. The split functions that occur are the same as for a simple 
impedance boundary condition, with each derivative factor in the boundary condition 
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giving rise to a pair of split functions. The edge diffraction coefficient then involves a 
product of these functions. 

The proposed boundary conditions have many possible practical 
applications. In analytical studies of the scattering from junctions and edges, they can 
be used to model single or multi-layered dielectric slabs or coatings. In numerical 
treatments, thick coatings can be simulated with a single condition on the surface of 
the layer, thus, reducing the required number of unknowns at the expense of a slight 
increase in the complexity of the integral equations. They could also be effective in 
modelling the dielectric layers used in printed circuit and microstrip arrays, leading to a 
significant simplification in the Green's functions involved, and we are now examining 
this possibility. 
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